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ABSTRACT 

We present the final band-merged ELAIS Catalogue at 6.7, 15, 90, and 175 

/im, and the associated data at U,g',r'-i'jZ,J,H.K, and 20cm. The origin 
of the survey, infrared and radio observations, data-reduction and optical 
identifications are briefly reviewed, and a summary of the area covered, 
and completeness limit for each infrared band is given. A detailed discus- 
sion of the band-merging and optical association strategy is given. The 
total Catalogue consists of 3762 sources. 23 % of the 15 /im sources and 
75 % of the 6.7 /im sources arc stars. For cxtragalactic sources observed 
in 3 or more infrared bands, colour-colour diagrams are presented and 
discussed in terms of the contributing infrared populations. Spectral en- 
ergy distributions (seds) are shown for selected sources and compared with 
cirrus, M82 and Arp220 starburst, and AGN dust torus models. 

Spectroscopic redshifts are tabulated, where available. For the Nl and 
N2 areas, the INT ugriz Wide Field Survey permits photometric redshifts 
to be estimated for galaxies and quasars. These agree well with the spectro- 
scopic redshifts, within the imcertainty of the photometric method (~ 10% 
in (1+z) for galaxies). The redshift distribution is given for selected ELAIS 
bands and colour-redshift diagrams are discussed. 

There is a high proportion of ultraluminous infrared galaxies {Igio of 
1-1000 /im luminosity L,r > 12.22) in the ELAIS Catalogue (14% of 15 
/im galaxies with known z), many with Arp220-like seds. 10 % of the 15 
fjja sources are gemiine optically blank fields to r' = 24: these must have 
very high infrared-to-optical ratios and probably have z > 0.6, so are high 
luminosity dusty starbursts or Type 2 AGN. 9 hyperluminous infrared 
galaxies {Lir > 13.22) and 9 EROs (r-K > 6) arc foimd in the survey. The 
latter are interpreted as ultraluminous dusty infrared galaxies at z ~ 1. 
The large numbers of ultraluminous galaxies imply very strong evolution 
in the star- formation rate between z = and 1. There is also a surprisingly 
large population of luminous {L^r > 11.5), cool (cirrus-type seds) galaxies, 
with Lir — Lopt > 0, implying Ay > 1. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: starburst - cosmology: observations 



The European Large Area ISO Survey (ELAIS) was 
originally proposed in response to the first Infrared 
Space Observatory {ISO) call for open time pro- 
posals in 1995 by M. Rowan- Robinson and 11 co- 
investigators from 9 institutions (Rowan- Robinson et 
al 1999). Subsequently the collaboration has grown 
to a total of 77 investigators from 32 European insti- 
tutions. The original concept was for a survey of 12 
sq deg of sky at wavelengths of 15 and 90 /tm. Sub- 
sequent awards of observing time allowed the survey 
to be extended to 6.7 /xm and (in collaboration with 
the FIRBACK team led by J.-L.Puget) 175 /im. The 
ELAIS areas were also surveyed at 20 cm with the 
VLA and the AT. 



1 INTRODUCTION 



which incorporates a physical models of cosmic ray 

and transient effects, was given by Lari et al (2001), 
and a first application of this to the ELAIS SI area 
was described by Gruppioni et al (2001). The final 15 
/tm reduction of the SI area is reported by Lari et al 
(2004), of the S2 area by Pozzi et al (2003) and of 
the Nl, N2 and N3 areas by Vaccari et al (2004) (see 
below for explanation of survey areas). The final re- 
duction of the 90 /im data is reported by Heraudeau et 
al (2004) and of the FIRBACK-ELAIS 175 /tm data 
by Dole et al (2002). A future paper by Rodighiero 
et al (2004, in preparation) will report results of ap- 
plying the Lari method to the ELAIS 90 /tm data. 
The reduction and analysis of the 20 cm data was de- 
scribed by Ciliegi et al (1999) for Nl, N2 and N3, and 
by Gruppioni et al (1999) for SI and S2. 



The survey goals, selection of survey areas, and 
details of the ISO observations, were described by 
Oliver et al (2000), and a preliminary analysis of the 
6.7 and 15 /tm data was reported by Serjeant et al 
(2000) and of the 90 /im data by Efstathiou et al 
(2000). A new method of reduction of the 15 /tm data, 



Associated with the ELAIS survey there has also 
been an extensive program of ground-based optical 
and near- infrared imaging and spectroscopy. The op- 
tical and spectroscopic follow-up of the SI area has 
been presented by La Pranca et al (2004), of the S2 
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area by Pozzi ct al (2003), and of the Nl and N2 ar- 
eas by Horaudeau at al (1999), Vaisanen et al (2002), 
Gonzalez- Solares et al (2004), Perez- Fournon et al 
(2004), Serjeant et al (2004), Verma et al (2004) and 
Afonso-Luis et al (2004). 

Deeper surveys at 6.7 and 15 /Ltm were also carried 
out by a subset of the ELAIS consortium in HDF-N 
(Serjeant et al 1997, Goldschmidt et al 1997, Oliver 
et at 1997, Mann et al 1997, Rowan- Robinson et al 
1997) and in HDF-S (Oliver et al 2002, Mann et al 
2002). 

X-ray surveys have also been carried out in sev- 
eral ELAIS areas. Alexander ct al (2001) reported 
Beppo-Sax observations over a large fraction of SI. 
Manners et al (2003) have reported Chandra obser- 
vations in the central regions of Nl and N2, and some 
interpretation of these observations has been given 
by Willott et al (2003). ROSAT data were compared 
with the ELAIS Preliminary Catalogue by Basilakos 
et al (2002). 

The present paper reviews the parameters of the 
ELAIS survey, gives a detailed account of the merg- 
ing of the different individual wavelength surveys to 
generate a final catalogue, and discusses the popula- 
tions present in the survey through colour-colour and 
colour-redshift diagrams, and spectral energy distri- 
butions. 



2 QUALITY OF THE CONSTITUENT 
INFRARED CATALOGUES 

The separate wavebands making up the ELAIS sur- 
vey each comprise an independent survey and are dis- 
cussed in separate papers. The present paper does not 
replace these analyses but pulls together those aspects 
which emerge from band-merging the surveys. For de- 
tailed discussion and analysis of the ELAIS sources, it 
is essential to refer back to these analyses of the indi- 
vidual constituent surveys. Table 1 defines the survey 
areas, wavelengths, characteristic depth, and source- 
densities. 

The calibration, completeness and reliability of 
the ELAIS 15 /im survey final analysis has been dis- 
cussed by Lari et al (2001, 2003), Gruppioni et al 
(2001) and Vaccari et al (2003). In the Lari method 
(Lari et al 2003), the ISOCAM Handbook (Blom- 
maert et al 1998) conversion factor from 750 instru- 
mental units (ADU) to mjy of 1.96 is used, with 
full modelling of the detector stabilization. To achieve 
consistency with the IRAS Faint Source Catalog 12 
/ttm calibration a multiplicative factor of 1.097 is ap- 
plied to 15 /X fluxes, based on the analysis by Aus- 
sel and Alexander (2003, in preparation) of near in- 
frared, ISO and IRAS fluxes for stars (Vaccari et al 
2003). The IRAS 12 //m calibration has been tested 
by Aussel and Alexander for over 4000 stars detected 
by IRAS and by 2MASS. By contrast, the Preliminary 
Analysis (PA) of Serjeant et al (2000) used a direct 



conversion of 1.75 ADU/mJy, based on a comparison 
of ISO and IRAS fluxes for stars, without correction 
for detector stabilization. Subsequently Vaisanen et 
al (2002) used J and K observations of ISO stars to 
derive a corrected conversion factor for the PA Cata- 
logue of 1.05 ADU/mJy. Figure 1 shows a plot of the 
final 15 iim flux, Sis, versus the PA flux (in ADU) 
(Babbedge 2004), which supports the conclusions of 
Vaisanen et al (2002). The Final Analysis Catalogue 
goes significantly deeper than the PA Catalogue, but 
almost all PA sources are found to be real, though 
much of the scatter in Fig 1 must be attributed to 
lower accuracy of the PA fiuxcs. 

It would also have been desirable to apply the 
Lari method to the 6.7 /im survey, but we have not 
been able to put together the very considerable re- 
sources needed to do this. As we will see in section 3.3 
75 % of the 6.7 /im sources are stars so the scientific 
returns for the present, predominantly extragalactic, 
study did not seem great. We have relied on the analy- 
sis of Vaisanen et al (2002) which confirmed both the 
reality of the 6.7 /xm Preliminary Analysis sources, 
and their flux calibration. The SIRTF-SWIRE survey 
(Lonsdale et al 2003) will allow comparison of the 6.7 
ixm ISO fluxes with SIRTF 8 /im fluxes, for sources in 
N2 and SI. The 6.7 /xm PA Catalogue has duplicate 
entries where different rasters observed the same part 
of the sky and only the entry with lowest noise was 
included. 

The calibration, completeness and reliability of 
the ELAIS 90 fim survey flnal analysis has been dis- 
cussed by Heraudeau et al (2004). Calibration of 90 
/xm photometry is based on a direct comparison of 
standard stars. Figure 2 shows a plot of the final 90 
/xm flux, 590, versus the PA flux (Babbedge 2004). 
Again, most of the PA sources arc confirmed. The 
substantial scatter can be mainly attributed to the 
limitations of the PA processing. 

The 175 /xm data analysis and identifications 
have been discussed by Dole et al (2002). A 5-cr sen- 
sitivity of 223 mJy was achieved. 

3 FINAL BAND-MERGED CATALOGUE 
3.1 Band-merging 

We have proceeded with the band-merging of the 
multi-wavelength ELAIS catalogues in a sequential 
way, taking into account the different positional ac- 
curacies of the component catalogues. The 1-ct po- 
sitional uncertainties at 15 /xm and 20cm have been 
estimated to be ~ 1" (Lari et al 2004, Vaccari et al 
2004, Gonzalez-Solares et al 2004, Ciliegi et al 1999, 
Gruppioni et al 1999), so sources at these wavelengths 
can be very reliably Eissociated with optical counter- 
parts, down to at least r' ~ 23 mag. We therefore first 
separately identify the 15 /xm and 20 cm catalogues 
with optical sources, using the likelihood method of 
Mann et al (1997), as discussed by Gonzalez-Solares 
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Table 1. Summajry of survey wavelengths, areas and numbers of sources, characteristic depth, source-densities 



Name 


RA 


Dec 


6.7 


15 


90 


175 f^ra 


Nl 


16''10'"01'' 


+54°30'36" 




2.67/490 


2.56/151 


2/103 (sq deg/no. of sources) 


N2 


16 36 58 


+41 15 43 


2.67/628 


2.67/566 


2.67/174 


1/55 


N3 


14 29 06 


+33 06 00 


1.32/189 


0.88/131 


1.76/119 




SI 


00 34 44 


-43 28 12 


1.76/304 


3.96/317 


3.96/226 




S2 


05 02 24 


-30 35 55 


0.12/40 


0.12/43 


0.11/5 






total area/number 




5.86/1161 


10.3/1546 


11.06/674 


3.0/158 sq deg/no. of sources 




characteristic depth 




1.0 


0.7 


70 


223 mjy 




source-density 




200 


150 


61 


53 per sq deg 




1 10 100 

15 micron prelim flux, S{PA), in ADU 



Figure 1. Comparison of Final Analysis 15 fim flux (mJy) 
(Lari et al 2004, Vaccari et al 2004) with Preliminary Anal- 
ysis 15 ixm fluxes (ADU) (Serjeant et al 2000). The Une 
corresponds to equal fluxes. 



S 

S 1000 



E 100 



100 1000 
90 micron prelim flux, S{PA), in mJy 



Figure 2. Comparison of Final Analysis 90 firn flux (mJy) 
(Heraudeau et al. 2003) with Preliminary Analysis 90 nm 
fluxes (Efstathiou et al 2000). 



et al (2004). The latter have described detailed sim- 
ulations of this association process, which show that 
the probability of spurious associations is < 5% for 
r' < 20, rising to 20 % by r' = 24. The 15 and 
20 cm sources are then merged on the basis of their 
optical positions, using a search radius of 2". If the 
optical positions agree to this accuracy, the radio flux 
and error are added to the 15 iim catalog entry. The 
percentages of 15 /um and 20 cm sources which found 
matches with the other wavelength were 8% and 11% 
respectively. A flag is set for sources for which there 
is more than one candidate identification (see below) : 



the selected association is the one with the highest 
likelihood (Gonzalez-Solares et al 2004). 

Radio sources which did not find a 15 iim match 
are then interleaved with the matched sources to give 
an RA-ordered list. Matched sources are given the 15 
/um source name. 

We next matched the 6.7 fim PA sources with 
the combined 15 /Lim-20 cm list, using a search radius 
of 5" . The larger search radius is required because of 
the poorer astrometry of the PA Catalogue, which is 
based on a single raster, with pixel size 3" . 30 % of 6.7 
/im sources found a match with a Catalogue source. 

Non-matched 6.7 /um sources were associated 
with optical counterparts where possible, and then 
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interleaved to generate a combined 6.7-15 pm-20 cm 
source-list. The positions of all sources in this list 
which did not have 20 cm fluxes were examined in 
the ELAIS 20 cm data to see whether a source in the 
3-5 (7 range might bo present, and also in the VLA 
FIRST survey for areas not surveyed by ELAIS, and 
if so the flux was added to the Catalogue. Otherwise 
3-cr radio limits are given as negative entries in the 
Catalogue. 

The 90 and 175 /im sources were then associ- 
ated with this combined list using search radii of 30 
and 60" respectively. Where a 90 or 175 /xm source 
is matched with more than one Catalogue entry, the 
less probable associations are flagged and are not used 
in subsequent discussions of infrared colours or spc- 
tral energy distributions, i.e. all the flux is assigned 
to the most likely match. The non-matched 5-(t 90 
and 175 /im sources are retained as a separate sup- 
plementary catalogue, because of their much poorer 
positions. They were searched for optical counterparts 
of high likelihood using search radii of 45 and 90" . Be- 
cause of the large positional uncertainties at 90 and 
175 ^ra, optical counterparts have to be reasonably 
bright to have a high likelihood of being the correct 
association (r < 19.0). 20 % of 90 /xm sources and 1 % 
of 175 n m sources, brighter than 5 a, found matches 
neither with the combined 6-7-15 /Lim-20 cm source- 
list nor with optical counterparts. Most of these may 
be associated with galaxies fainter than r' = 19. How- 
ever we can not rule out the possibility that some of 
the unassociated 90 yum sources may be spurious. 

Finally the combined 6.7-15-90-175 /im-20 cm 
Catalogue is searched for pairs within 5" and each 
of these is examined individually. In most, but not 
all, cases these pairs are believed to be the same ob- 
ject, split in two by the optical association process, 
and in these case the sources have been merged, with 
an appropriate flag set in the Catalogue. 

3.2 Catalogue Description 

The Final Band-merged Catalogue is given at 
http://astro.imperial.ac.uk/Elais/ . The Catalogue 
entries are : source name (for merged sources, in or- 
der of preference 15 pm, 20 cm, 6.7, 90, 175 fim. ), 
source position (same order of preference for merged 
sources), 20 cm flux and error, 175 jjxa flux, error, 
S/N, positional ofi^set, 90 /um flux, error, S/N, posi- 
tional offset, 15 /im flux and S/N, 6.7 /im flux and 
error, 6.7 /xm reliability flag (2 = high reliability, 3 = 
medium reliability), flag 1 (see below), flag 2, J mag 
and error, H mag and error, K mag and error, flag 
3, position of optical association, WFS U, g', r', i', 
Z magnitudes and errors, WFS star/galaxy flags (-1 
for stellar image, 1 for extended image), Sextractor 
r' mag and error, Sextractor star/galaxy classificar 
tion (0.0 for galaxy, > 0.7 for stellar image), posi- 
tional offset of optical ID (total, RA, Dec), probabil- 
ity of optical association (threshold 0.8, see eqn (5) 



of Gonzalez-Solarcs ct al (2004) for definition), relia- 
bility of optical association (see eqn (6) of Gonzalez- 
Solares et al (2004) for definition), photometric spec- 
tral energy distribution (sod) type, ntyp, and redshift 
(lgio(l + Zphot), derived assuming Av = 0), photo- 
metric sed type and redshift and Ay (free fit for Ay), 
spectroscopic redshift Zspect, flag 4, flag 5, Zbest = 
Zspect if available, = Zphot otherwise, bolometric opti- 
cal luminosity Lopt, ir sed type (1 = cirr, 2 = M82-sb, 
3 = A220-sb, 4 = AGN dust torus, 5 = cooler cirrus, 6 
= 2-1-4), bolometric (1-1000 /Ltm) infrared luminosity 
Lir, AGN dust torus luminosity Ltor (if 15 /xm emis- 
sion is interpreted as due to dust torus emission). 

where flag 1 = 1 if radio flux force-merged (15 
um and radio positions within 5"), = 2 if 6.7 fj,m flux 
force-merged, = 3 if 90 fj,m association is not the most 
likely, where there are multiple associations, = 4 if 175 

fim association is not the most likely, where there are 
multiple associations, = 5 if second most probable 
optical association for 15 um has been preferred on 
basis of radio position, = 6 if source has an IRAS 
association, = 7 if 1 and 6 both set, = 8 if 3 and 4 
both set, = 9 if 1 and 8 both set, or 4 and 5 both set; 

flag 2 = 1 if source falls in gaps between WFS 
chips, so no optical data, = 2 if source falls near edge 
of WFS chip, so photometry may be unreliable, = 
3 if source hcis multiple optical counterparts, = 4 if 
source is blank in optical (ie no optical counterpart 
with probability of association > 0.7, within specified 
search radius), = 6 if there is an association in NED 
(non-ELAIS, non-2MASS, non-IRAS), = 7 if associa- 
tion is bright star, = 8 if flags 1 and 7 set, = 9 if flags 
3 and 7 set; 

flag 3 = 1 if J,H,K magnitudes from 2MASS, = 

2 if J, K magnitudes from Vaisanen et al (2002), = 

3 if K magnitudes from Rigopoulou et al (2004, in 

preparation), = 4 if K magnitude from Pozzi ct al 
(2003a) = 5 if K magnitude from Sajina et al (2003); 

flag 4, Uztyp (spectroscopic classification), = 1 
(spiral) galaxy, — 2 emission-line, starburst, — 3 ab- 
sorption line, early type, = 4 AGN, — 5 Syl, = 6 Sy2, 
= 7 star, =8 liner; 

flag 5, 1 Perez-Fournon et al (2003), = 

2 Serjeant et al (2003), = 3 La Franca et al (2003), 
= 4 Pozzi et al (2003a), = 5 NED, = 6 Morel et al 
(2001), = 7 SLOAN Survey, = 8 Willott et al (2003), 
= 9 Chapman et al (2002), Sajina et al (2003). 

For extended objects in Nl, N2, the WFS mag- 
nitudes refer to the flux within a fixed aperture: 
they should give the correct colours of the objects, 
as required by the photometric redshift code. For 
integrated magnitudes the Sextractor r' magnitude 
should be used (and other WFS bands can be cor- 
rected by the difference between the WFS and Sex- 
tractor r' magnitudes). Because of the effects of sat- 
uration, colours can only be trusted if the Sextractor 
r' magnitude is > 15 (this mainly affects stars in the 
Catalogue). In SI, the B magnitudes are derived from 
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APM magnitudes and are very uncertain (also possi- 
bly too faint on average by ~ 0.5 mag.). 

There are 3523 sources in the 6.7-15/im-20cm 
Catalogue, of which 1636 are 15 /xm sources, 1136 
are 20 cm (non 15 ^m) sources, 741 are 6.7 fim (non 
15 iim, non 20 cm) sources, and 239 sources are in the 
supplementary 90-175/im catalogue. The numbers of 
15 pm sources with redshifts are 204 in SI (199 spec- 
troscopic), 31 in S2 (22 spectroscopic), 10 in N3 (all 
spectroscopic), 309 in Nl (109 spectroscopic) and 355 
in N2 (167 spectroscopic), 909 15 /jm redshifts in all, 
and 1210 redshifts in the whole Catalogue (see below 
for discussion of photometric redshifts). 97 % of 15 
^m sources in Nl and N2 arc accounted for as cither 
stars, extragalactic sources with redshifts, or blank 
fields (8 and 11 % respectively), and 93 % in SI (the 
percentage of blank fields in SI is higher, 20 %, be- 
cause of the shallower optical survey). 48 % of the 
radio (non 15 jj,m) sources are blank fields (to r' = 
24). 91 % of the 6.7 fim (non 15 /xm or radio) sources 
in Nl are found to be stars. 

3.3 Associations 

We matched the whole Catalogue with the NASA 
Extragalactic Database (NED) and redshifts result- 
ing from these associations are included in our Cata- 
logue. We also specifically matched the Catalogue to 
the 2MASS J,H,K catalogues, finding matches for 30 
% of our sources (with a higher success rate for the 
stars). A search radius of 5" was used, and in the few 
cases where sources appeared in both the 2MASS Ex- 
tended and Point-Source Catalogues, the magnitudes 
from the Extended Catalogue were preferred. Here we 
discuss some more specific results from associations 
with known objects. 

Table 2 hsts the 750 and IRAS data for ELAIS 
sources which are associated with IRAS Faint Source 
Catalog sources (FSC, Moshir ot al. 1992), with up- 
per limits indicated as negative values. Of the 39 as- 
sociations in Nl and N2, 12 are clearly stellar photo- 
spheres detected by IRAS at 12 fim (and occasionally 
at 25 /um) and the rest are nearby normal galaxies. For 
the galaxies detected by IRAS at 60 and 100 fim, we 
have interpolated to estimate a 90 pm flux and calcu- 
lated logioS{90)iRAs/S{90)iso- The mean value for 
19 sources is 0.129 ±0.022, suggesting that there are 
unresolved calibration issues between IRAS and ISO. 
However it should be noted that there is a tendency 
for the IRAS FSC to overestimate fluxes near the FSC 
threshold at 60 and 100 /^m (Moshir et al 1992). Since 
all the IRAS sources in Table 2 have S(IOO) < 3 Jy 
(and most have S(IOO) < 1 Jy), this would be suffi- 
cient to explain the discrepancy noted above. 

The ELAIS Nl field was also partially observed 
with the Ha survey of Pascual et al. (2001). Since 
the infrared and Ha both trace star-formation it was 
tempting to see if there were any sources in com- 
mon. Using the Ha, FIR, star-formation calibrations 



of Cram et al. (1998) we were able to estimate a mean 
Ha flux to FIR flux. Then using a starburst SED 
we were able to estimate the expected 15/im flux for 
each of the Pascual et al. sources. With one exception 
the expected 15/im fluxes all fell below the charac- 
teristic depth of the survey (0.7 mjy. Table 1). The 
exception was the source in their field a3 with ID 
7227b at 16 05 46.3 -|-54 39 11.74 with mi = 17.5 
and mua = 17.0. We estimate that this source should 
have had a 15/ini flux of 2.9 m,]y. Neither this source 
nor any of the others were detected in our ELAIS 
Catalogue. A more detailed analysis of the expected 
dispersion in the Ha/FIK relation is required before 
we can assess whether these non-detections suggest 
that the mean relation needs to be revised. 

The Canada-France-Hawaii Telescope blue grens 
quasar survey of Crampton et al. (1992) overlaps with 
the N2 field. 79 of their candidates fall within the 
ELAIS boundaries. 11 of these candidates are de- 
tected by ELAIS and are listed in Table 3. 

Stars can be recognized in the Catalogue through 
(i) having flag 2 set to 7, (ii) having flag 4 set to 7, 
(iii) having low ratios of {Si^/Sr) or (Sa.r/Sr), where 
Sr is the r-band flux in mJy. Almost all sources with 
logioiSis/Sr) < —0.5, or logio{Sfi.7 / Sr) < —0.2, are 
stars. A total of 846 stars have been identified in this 
way in the Catalogue, comprising 22, 23 and 25 % 
of the 15 fim sources in N2, Nl and SI respectively, 
and 75 % of 6.7 /jm sources. In Nl, 91 % of sources 
which are 6.7 /xm (non-15 /xm or 20 cm) sources are 
stars. The infrared fluxes are consistent with being 
photospheric emission in almost all cases. The small 
number of stars where the ISO emission appears to be 
in excess of the photspheric prediction deserve further 
detailed study, but this is beyond the scope of the 
present paper. 



4 RADIO-INFRARED-OPTICAL 
COLOUR-COLOUR DIAGRAMS 

For sources detected in 3 bands we can plot colour- 
colour diagrams and compare the results with pre- 
dictions of models. Figs 3-8 show a selection of these, 
with predicted loci corresponding to the basic infrared 
templates used by Rowan-Robinson (2001): cirrus, 
M82 starburst, Arp220 starburst, AGN dust torus. 

Figure 3 shows the 175/90/15 colour-colour di- 
agram for Nl and N2, with loci for cirrus (C), M82 
starburst (S), and Arp220 starburst (A), for redshifts 
from to 3 (the position of the labels indicates the 
zero redshift end of the loci) . The cross in the lower 
left-hand corner of this and subsequent diagrams in- 
dicates typical (median) error bars (where these are 
larger than the plotted symbols). The model loci loop 
around the diagram so without additional redshift in- 
formation it is hard to say much about the popu- 
lations (or derive photometric redshifts from far in- 
frared data). However from both spectroscopic and 
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Table 2. Matches with the IRAS Faint Source Catalog Version 2 



IRAS name 


12 


25 


60 


100 fj,m 


IRAb-lbO 


ELAIS name 


6.7 


15 


90 


175 /J 




mjy 








sepn (') ) 


mJy 










F00279-4253 


-75.4 


-76.4 


649.0 


1440.0 


0.70 


ELAISC15-J003022-423657 


- 


25.21 


849 


- 


F00302-4249 


-112.0 


-99.0 


232.0 


-525.0 


0.53 


ELAISC15-J003244^423313 


- 


11.64 


191 


- 


F00315-4421 


178.0 


-93.8 


-123.0 


-376.0 


0.43 


ELAISC15-J003402-440442 


- 


109.34 


- 


- 


r 00320-4342 


-65.8 


-103.0 


309.0 


685.0 


0.25 


ELAloC15-J003429-432614 


6.18 


24.32 


408 




F00325-4313 


-96.6 


-75.5 


322.0 


906.0 


0.13 


ELAISC15-J003458-425733 




15.83 


534 




F00341-4428 


-129.0 


-155.0 


202.0 


-676.0 


0.78 


ELAISC15-,1003626-'111140 


5.80 


13.61 


290 


_ 


F00353-4418 


-53.4 


-70.8 


204.0 


-559.0 


0.20 


ELAISC15-J003741-440226 


_ 


1.36 


210 


_ 


F00360-4355 


-124.0 


-128.0 


374.0 


1280.0 


0.27 


ELAISC15-J003828-433848 


10.14 


46.58 


682 


_ 


F00362-4416 


135.0 


-149.0 


-166.0 


-279.0 


0.12 


ELAISC15-J003836-440029 


- 


27.48 


- 


- 


F14262-I-3328 


-60.2 


-59.7 


317.0 


518.0 


0.27 


ELAISC15-J142823.4+331513 


- 


12.72 


448 


- 


F14266+3336 


-67.7 


-98.9 


452.0 


999.0 


0.32 


ELAISC15- J142847. 1+332315 


_ 


19.02 


540 


_ 


F14292+3318 


143.0 


-116.0 


-105.0 


-316.0 


0.18 


ELAISC15-J143123.5+330517 


258.75 


62.91 


- 


- 


b 14292+3327 


-68.4 


94.9 


694.0 


983.0 


0.13 


ELAloU15- J 143125.3+331348 


6.93 


28.95 


933 




F14304+3341 


-68.1 


-47.6 


201.0 


476.0 


0.13 


ELAISC15-J143234.9+332833 




1.57 


434 




F16022+5450 


-45.0 


-65.7 


159.0 


-614.0 


0.70 


ELAISC15-J160322.8+544237 


- 


2.36 


101 


- 


F16029+5506 


-66.7 


-67.0 


296.0 


507.0 


0.09 


ELAISC15- .1160408.4+545812 


- 


8.77 


315 


- 


F16046+5415 


-83.2 


78.8 


604.0 


798.0 


0.25 


ELAISC15-J160552. 5+540650 


- 


30.09 


756 


838 


F16063+5405 


-69.3 


-57.6 


244.0 


-742.0 


0.16 


ELAISC15-J160736.5+535731 


- 


19.15 


230 


597 


F16070+5439 


97.5 


-50.4 


-93.2 


-308.0 


0.12 


ELAISC15-J160812.7+543141 


- 


53.25 


- 


- 


F16083+5400 


-48.3 


-42.6 


184.0 


-720.0 


0.11 


ELAISC15-J160934.7+535220 


_ 


1.99 


215 


309 


F16091+5357 


207.0 


69.9 


-102.0 


-419.0 


0.09 


ELAISC15-J161019.3+534934 


- 


103.67 


- 


- 


F16091+5447 


576.0 


146.0 


-91.2 7 


-244.0 


0.02 


ELAISC15-J161017.6+543929 


- 


288.39 


- 


- 


F16145+5447 


-35.7 


-64.4 


161.0 


475.0 


0.33 


ELAISC15-J161545. 8+544019 


- 


20.70 


303 


- 


F16294+4115 


-87.8 


-65.4 


243.0 


936.0 


0.13 


ELAISC7-J163104+410913 


2.96 


- 


540 


- 


F16298+4129 


131.0 


-53.5 


-124.0 


-424.0 


0.22 


ELAISC15-J163130.2+412330 


255.16 


66.21 


- 


- 


F16323+4127 


-74.9 


-67.8 


383.0 


743.0 


0.23 


ELAISC15-J163401.8+412052 


9.60 


20.37 


403 


666 


F16334+4116 


-48.4 


-65.0 


190.0 


-813.0 


0.17 


ELAISC15-J163506.1+411038 


- 


7.98 


251 


346 


F16337+4101 


-66.4 


-73.7 


224.0 


-843.0 


0.17 


ELAISC15-J163525. 1+405542 


5.46 


13.91 


416 


682 


F16338+4138 


147.0 


-70.2 


-104.0 


-322.0 


0.12 


ELAISC15-J163531. 1+413158 


524.98 


106.80 






F16341+4053 


103.0 


-80.2 


-88.2 


-369.0 


0.24 


ELAISC15-J163549.0+404720 


196.97 


58.42 






F16341+4059 


243.0 


-109.0 


-102.0 


-339.0 


0.10 


ELAISC15-J163549.9+405317 




125.86 






F16344+4111 


-71.9 


-67.6 


351.0 


897.0 


0.06 


ELAISC15-J163608.1+410507 


2.30 


8.94 


614 


803 


F16349+4038 


-84.7 


-35.8 


2.36.0 


615.0 


0.37 


ELAISC15-J163633.5+403245 


6.15 


16.52 


291 




F16349+4034 


116.0 


-50.7 


-98.7 


-441.0 


0.04 


ELAISC15- .1163637.3+402824 




49.44 






F 16359+4058 


-61.1 


104.0 


1220.0 


2480.0 


0.19 


ELAISC15-J163734.4+405208 


18.37 


53.56 


1461 


2377 


F16365+4202 


-69.6 


84.4 


460.0 


1500.0 


0.28 


ELAISC15-J163814.0+415620 


6.14 


50.15 


611 




F16377+4150 


-54.3 


-68.8 


306.0 


520.0 


0.16 


ELAISC15-J163924.0+414442 


2.55 


8.13 


299 




F16389+4146 


249.0 


81.2 


-86.2 


-500.0 


0.06 


ELAISC15-J164033.9+414107 


645.30 


125.86 






F16405+4113 


-55.6 


-68.4 


169.0 


-545.0 


0.15 


ELAISC15-J164211.9+410816 


1.80 


8.32 


288 





photometric redshifts (see section 5) we can deduce 
that almost all these galaxies have z < 0.5. Most of 
the objects then lie reasonably close to the cirrus lo- 
cus over this redshift range, so we deduce that this 
bright subset consists of normal spirals with far in- 
frared emission from interstellar dust bathed in the 
general stellar radiation field. 

Figure 4 shows the 90/15/6.7 colour-colour dia- 
gram for N2 and SI, with the same model loci. Again 
there is ambiguity about the populations involved, 
which is not resolved using the fact that most galax- 
ies detected at 90 jim have z < 0.5. An earlier version 
of this figure was discussed by Marquez et al (2002). 

Figure 5 shows the 90/15/r-band colour-colour 
diagram for Nl, with the same model loci. Here the 



model loci are more differentiated and we can deduce 
that all three templates are represented, with cirrus 
and M82 starburst components accounting for over 75 
% of the sources. 

Optical colours provide a powerful discriminant 
between galaxies and AGN. Many of the Catalogue 
sources in Nl and N2 are detected in each of the WFS 
band g', r', i'. In figure 6 we show the g' —r' —i' colour- 
colour diagram for sources classified by Sextractor as 
galaxies. The brighter galaxies define a very tight set 
of colours, narrower in (r' — i') than in (g'-r)', char- 
acteristic of galaxies with < z < 0.5. Model curves 
are shown for E, Sbc and starburst galaxies, for z = 
0-2. In figure 7 we show the corresponding plot for 
objects classified as star-like (excluding actual stars), 
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Table 3. Quasar candidates from the blue grens survey of Crampton et al. 1992 that are detected by ELAIS. The name of 
the sources in the Crampton et al catalogue and the ELAIS Catalogue are given. This is followed by the separation of the 
two in arc seconds. The fourth column gives the spectral types estimated from the U, g,r,i,Z INT data in our photometric 
redshift fitting. The Q flag indicates the confidence of the quasar candidate in the Crampton et al. work from 1 (strong 
lines) through to 4 (blue continuous spectrum). 



Crampton et al. name 


ELAIS name 


Sep/arc sec opt s« 


!d type 


Q 


1635.20+4124 


ELAISC15 J163652.7+411827 


1.2 


4 


2 


1635.30+4135 


ELAISC15 ,1163659.0+412928 


2.5 


6 


3 


1635.40+4136 


ELAISC15 ,1163702.2+41.3022 


1.1 


7 


4 


1635.70+4121 


ELAISC15 J163721. 3+411503 


0.7 


8 


3 


1636.00+4149 


ELAISC15 J163739.3+414348 


1.6 


7 


2 


1636.50+4203 


ELA1SC15 J163805.6+415740 


1.4 


3 


3 


1636.70+4133 


ELAISR163817+412730 


1.3 


7 


2 


1637.20+4217 


ELAISC15 J163847.5+421141 


0.6 


7 


4 


1637.60+4134 


ELAISC15 J163915.9+412834 


1.0 


6 


3 


1638.60+4126 


ELAISC15 J164016.0+412102 


2.5 


7 


1 


1638.70+4108 


ELAISC15 J164018.8+410254 


1.2 


8 


3 
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Figure 3. 175-90-15 fj,m colour-colour diagram for sources 

in Nl and N2, with loci for cirrus (C), starburst (S) and 
Arp220 (A) scds, from z =0 to 3 (labels denote the z = 
end of the locus, crosses mark z = 1, 2, 3). Only sources 
detected in all 3 bands are plotted: all are galaxies. The 
cross in the lower left-hand corner of this and subsequent 
figures indicates typical (median) error bars. 

with different symbols for sources which the photo- 
metric redshift code (sec section 5) classifies as having 
AGN-typc seds and those with galaxy seds. The AGN 
occupy quite a tight colour region centred on (g — r) ^ 
0.4, (r' — i') ~ 0.3. Model curves are shown for E and 
starburst galaxies for z = 0-2 and for AGN with z 
= 0-6. Objects with galaxy seds show some overlap 
with Fig 6, but with more scatter to higher values of 
(r' — i'), consistent with having higher redshifts. 

Figure 8 shows lg(S'i5/5r) versus (w — r) for 
ELAIS sources , where Sr denotes the r'-band flux in 
mjy. There is a clear separation between the Galactic 



Figure 4. 90-15-6.7 /im colour-colour diagram for sources 

in N2 and 81 detected in all 3 bands (all are galaxies), with 
loci for cirrus (C), starburst (S) and Arp220 (A) seds, from 
2 = to 3. 

stars in the lower part of the diagram and the AGN 
and compact galaxies in the upper paxt, and between 
the galaxies and AGN. The objects classified as star- 
like but with galaxy seds have values of Si^/Sr up 
to 300 and these must be heavily obscured starbursts 
hke Arp 220, or Type 2 AGN. 



5 PHOTOMETRIC REDSHIFTS IN Nl 
AND N2 

The WFS optical data in the U,g',r',i',Z bands, and 
J,H,K data (where available) for Nl and N2 allows us 
to determine photometric redshifts for a large fraction 
of the sample. The approach used is that of Rowan- 
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Figure 5. 90-15-r ^vci colour-colour diagram for Nl , N2 
and SI, with loci for cirrus (C), starburst (S) and Arp220 
(A) seds, from z = to 3. 



Robinson (2003) with a set of 6 galaxy templates, and 
with the option of varying Ay- In addition two sim- 
ple AGN templates are included, based on the aver- 
age optical QSO spectrum of Rowan- Robinson (1995), 
modified to take account of observed seds of ELAIS 
AGN (details of all the templates used are given at 
http://astro.ic.ac.uk/~mrr/photz). Figs 9, 10 show 
the comparison of the data for AGN with spectro- 
scopic redshifts with the assumed templates. The ap- 
plication of this code to the full WFS data set is do- 
scribed by Babbedge et al (2004). We have also used 
the UBRI data in S2 (Pozzi et al 2003) to estimate 
photometric redshifts. In SI wc have only BRI data 
in the optical and the photometric calibration at B 
is uncertain. We found that without U (or u) data, it 
was impossible to determine photometric redshifts for 
quasars, but results for galaxies were still good, and 
these are included in the Catalogue provided at least 
3 photometric bands between B and K are available. 
Figure 11 shows the comparison of photometric and 
spectroscopic redshifts for galaxies in Nl, N2, SI and 
S2. The spectroscopic redshifts in Nl and N2 are re- 
ported by Perez-Fournon ct al (2004) and Serjeant et 
al (2004), in SI by La Franca et al (2004) and in S2 
by Pozzi et al (2003). The agremeent is good, within 
the uncertainty of the photometric method (~ 10% in 
(1-l-z) according to the analysis of Rowan-Robinson 
(2003)). A fuller discussion of the application of pho- 
tometric redshift techniques to the WFS data is given 
by Babbedge et al (2004). Of course, it needs to be 
emphasized, that a 10 % accuracy in (l-|-z) means 
that photometric redshifts << 1 will be very inaccu- 
rate. One of the two strongly discrepant points in Fig 



Figure 6. g' — t' — i' colour-colour diagram for galaxies 
in Nl and N2. Filled circles: r' < 21 mag., open circles: r' 
> 21 mag. Model loci for E (solid line), Sbc (dotted line) 
and starburst (broken line) galaxies are shown for z = 0-2 
(crosses denote z = 0.5, 1, 1.5, 2). 



1.5 



0.5 - 



- 



-0.5 



1/ 



g-r 



Figure 7. g' — r' — i' colour-colour diagram for star-like 
objects in Nl and N2. Filled circles: AGN sed, open circles: 
galaxy sed. Model loci for a quasar (long-dashed line, z = 
0-6), E galaxy (solid line, z = 0-2), starburst (dashed line, 
z = 0-2), are shown (crosses: z = 0.5, 1, 1.5, 2). 
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Figure 8. 15/r versus {U — r') for ELAIS galaxies and 
AGN in Nl and N2. Filled circles: extended optical identi- 
fications, filled triangles: star-like objects with AGN seds, 
open circles; star-like objects with galaxy seds, crosses: 
Galactic stars. Note that the (U-r') colours of the Galac- 
tic stars are inaccurate because of saturation effects. 



11 is a case where the photometric redshift is deter- 
mined from only 3 photometric bands, so as in Rowan- 
Robinson (2003) we can claim excellent performance 
for the photometric redshift method if 4 bands are 
available. Figure 12 shows the corresponding plot for 
AGN in Nl, N2 and S2: the results are surprisingly 
good. Several AGN acquire spurious photometric red- 
shifts in the range 2-2.5 because their U magnitudes 
are fainter than predicted by the template, possibly 
because of the effect of extinction, so they are inter- 
preted as Lyman drop-outs. The code is quite suc- 
cessful both in finding the AGN and estimating their 
redshifts. Out of 33 sources with both spectroscopic 
and photometric redshifts, which are classified either 
spectroscopically or photometrically as AGN, only 2 
are not classified spectroscopically as AGN, presum- 
ably due to aliassiug in the photometric solution. 5 
spectroscopic AGN are not picked up as AGN pho- 
tometrically though all have consistent photometric 
redshift estimates, perhaps due a weak AGN contin- 
uum. Examples of the template fits to optical and 
near ir data can be seen in Figs 26-28. 

Figures 13-16 show redshift histograms at 15, 90 
and 175 /um and 20 cm, for sources brighter than 
the characteristic depth specified in Table 1, with 
the dotted lines in Figs 11 and 14 indicating spectro- 
scopic redshifts. The broken histogram indicates the 
effect of assigning the 15 /xm blank-field sources in SI 
(which have R > 20) redshifts uniformly distributed 
in the range 0.2-1.5, and the blank-field sources in 



Figure 9. Comparison of photometric data for AGNs with 
ntyp = 7 and known spectroscopic redshifts with the as- 
sumed template. 




Figure 10. Comparison of photometric data for AGNs 
with ntyp = 8 and known spectroscopic redshifts with the 
assumed template. 



Nl and N2 (which have r' > 24) redshifts uniformly 
distributed in the range 0.6-1.5, based on the R-z di- 
agram of LaFranca et al (2003, their Fig 15). The 
median redshift is 0.30 at 20 cm, 0.17 at 15 fxm and 
0.10 at 90 and 175 /xm. In their analysis of the 15 /xm 
redshift distribution in SI, La Franca et al (2003) and 
Pozzi et al (2004) infer a stronger secondary peak at 
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Figure 11. lgio(l+%hot) versus lgm(l+Zspect) for galax- 
ies. The straight Hnes bracket the range Algio(l + z) = 
±0.1 (±2.5f7 according to the analysis of Rowan-Robinson 
(2003)). The photometric redshifts are The straight lines 
indicate the ±2.5<t range. The photometric redshifts are 
for an assumed Ay = 0. 



Figure 13. Rcdshift histogram for ELAIS 15 ^nxi sources. 
Solid curve: both photometric and spectroscopic redshifts 
(840 sources), dotted curve: spectroscopic redshifts only 
(468 sources), broken curve: effect of assigning blank fields 
uniformly to range 0.5 < z < 1.5. Long broken curve: pre- 
dicted redshift distribution from Rowan-Robinson (2001), 
dash-dotted curve: predicted redshift distribution from 
Pozzi et al (2004). 




0.2 0.4 0.6 0.8 

lg(l- + Zsp.ct) 



Figure 12. Ig^oi^+^phot) versus lgio(l+Zspect) for AGN. 



around z = 1 than is indicated here. Model predic- 
tions by Rowan-Robinson (2001) at 15, 90 and 175 
/um are shown in Figs 13-15, and by Pozzi et al (2003) 
at 15 /xm in Fig 13. The overall agreement with the 
Rowan-Robinson (2001) model is reasonable, though 



the Pozzi et al (2003b) model may be a slightly better 
fit at 15 /xm. This merits further investigation. 

The luminosity function and evolution at 15 /xm 
are discussed by La Franca et al (2004), Pozzi et al 
(2004), Percz-Fournon ct al (2004), at 90 /xm by Ser- 
jeant et al (2004) and at 175 /jm by Perez-Fournon et 
al (2004). 

Figure 17 shows lg(S'i5/5'r) versus z, where z = 
Zbest, for objects classified by Scxtractor as galax- 
ies, together with predicted loci for cirrus, M82 and 
Arp220 starbursts. Although most of the galaxies 
have z < 0.3, there are an interesting subset with 
0.6 < z < 1.2. Figure 18 shows the corresponding di- 
agram for star-like objects, with model loci for AGN 
dust tori, and M82 and Arp220 starbursts. The ob- 
jects with AGN seds, shown as filled circles, follow 
the (Type 1) AGN dust torus model line well. Most 
with galaxy seds (and some from Fig 17 with 0.2 < 
z < 0.8) have Si^/Sr values higher than any of the 
model loci. These may represent a new population of 
heavily obscured starbursts or Type 2 AGN. 

Fig 19 shows \g{Sir, / S2()cm) versus z for all 
sources, for an assumed rest-frame Seo/^ocm ratio 
of 200. There is a dependence of Si5/S20cm on in- 
frared sed type and on redshift. Thus even if there is 
a perfect correlation between 60 jim and 20 cm fiux, 
we do not expect this correlation to be preserved at 
other infrared wavelengths. Note the two radio-loud 
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Figure 14. Redshift histogram for ELAIS 90 sources. 
Solid curve: both photometric and spectroscopic rcdshifts 
(229 sources), dotted curve, spectroscopic redshifts only 
(175 sources), Long broken curve: predicted redshift dis- 
tribution from Rowan-Robinson (2001). 



AGN in the lower part of Fig 19. The 15 pm - radio 
correlation has also been discussed by Gruppioni et 
al (2003). 

Fig 20 shows lg(S'i5/S'6.7) versus redshift for all 
sources in N2 and SI. Values of Sn/Se.v < 1 are pre- 
dominantly stars and have been correctly recognized 
as such by the photometric redshift code. The 15/6.7 
colour ratio is not a good discriminant between the 
gala^xy models at low redshift, though it is clear that 
at high redshift the sources conform well to the AGN 
dust torus locus. 

Figures 21 and 22 show lg(5'9o/S'r) and 
lg('S'i75/5'r) versus z, with model loci for cirrus, M82 
and Arp220 starbursts. The models nicely bracket the 
observational points. However a plot of lg(S'i75/S'go) 
versus lg(l + z) (Fig 21) shows that at low redshifts 
there is a group of nearby galaxies with colder dust 
( vBu{T) colour temperatures in the range 18-25 K) 
than the cirrus template used here (which matches 
IRAS colours for cirrus sources well). Clearly we do 
not expect all the dust in a galaxy to be at a sin- 
gle temperature. As we move to longer wavelengths 
we expect to see emission from cooler dust further 
from the centre of the galaxy. Observations of these 
galaxies at 350-850 /xm would be valuable for under- 
standing the true distribution of dust mass and tem- 
perature in galaxies. 

Figure 24 shows (w — r) versus lgio(l + «) for 
objects classified as extended, compared with model 
predictions for E, Sbc and starburst galaxy seds. The 
photometric redshift code suggests that galaxies with 



Figure 15. Redshift histogram for ELAIS 175 ^m sources. 
Solid curve; both photometric and spectroscopic redshifts 
(46 sources). Long broken curve; predicted redshift distri- 
bution from Rowan- Robinson (2001). 

(u-r) > 2 show evidence of extinction, with Av in the 
range 1-3. Figure 25 shows the corresponding diagram 
for star-like objects. The objects with AGN-type seds 
follow a narrow locus, with objects with spectroscopic 
redshifts (filled circles) agreeing well with the adopted 
AGN scd, while those with galaxy seds occupy a sim- 
ilax region to the galaxies of Fig 24. 

A full discussion of the JHK, 6.7 and 15 fjcm 
colour-colour diagrams has been given by Vaisanen et 
al (2002) and the results arc not significantly altered 
by the larger and more accurately calibrated sample 
here. The JHK data help the photometric redshift 
fitting but do not help understand the nature of the 
dust-emission components in galaxies. 



6 INFRARED SPECTRAL ENERGY 
DISTRIBUTIONS 

For sources detected in multiple ISO bands we can 
compare the infrared spectral energy distributions 
with those for standard model templates. Figure 26 
compares the rest-frame seds of 9 galaxies detected 
in all four ISO bands with a standard cirrus emis- 
sion spectrum (Efstathiou and Rowan- Robinson 2003: 
surface brightness parameter = 5, age of star- 
burst = 5 Gyrs). All have modest infrared luminosi- 
ties {Lir < 11.5). Figure 27 shows a similar plot for 
luminous cirrus galaxies {Lir > 11.5). These repre- 
sent an interesting new population of luminous cool 
galaxies, with redshifts in the range 0.15-0.5, con- 
sistent with the strong evolution postulated for cir- 
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Figure 16. Rcdshift histogram for ELAIS radio sources. 
Solid curve: both photometric and spectroscopic redshifts 
(424 sources), dotted curve: spectroscopic redshifts only 
(140 sources). 



Figure 18. 15/im/r-band colour versus z for ELAIS star- 
like sources (excluding stars), with loci for AGN (Q, solid 
line), starburst (B, dotted line) and Arp220 (A, broken 
line) seds. Filled circles: sources with AGN optical seds, 
open circles: sources the galaxy seds. 
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Figure 17. 15^m/r-band colour versus z for ELAIS galax- 
ies , with loci for cirrus (C, solid line), starburst (B, dotted 
line) and Arp220 (A, broken line) seds. 



rus galaxies by Rowan-Robinson (2001) and perhaps 
related to the very luminous cirrus galaxies postu- 
lated by Efstathiou and Rowan-Robinson (2003) to 
explain some of the high redshift SCUBA galaxies. 
For 160552.5-1-540650 there is evidence also for a star- 
burst component. Figure 28 compares 2 galaxies de- 



tected in 4 ISO bands and 4 detected in 3 ISO bands 
with a high optical depth starburst model which gives 
a good fit to Arp 220 (Efstathiou et al 2000: Ay = 
200, age of starburst = 26 Myr). The tend to have 
higher infrared luminosities and optical sed type 
= 1 or 2 (E or Sab), consistent with the optical sed 
of Arp220, which shows very little contribution from 
young massive stars duo to the high extinction in this 
gala^xy. Figure 29 compares 3 gala^Kies detected in 3 
ISO bands with an M82-type starburst model (Efs- 
tathiou and Rowan-Robinson 2001: Av = 50, age of 
starburst — 26 Myr). Two galaxies for which an AGN 
dust torus component is also required are included, 
as well as a further two objects which only require an 
AGN dust torus component. These figures include all 
ultraluminous ELAIS galaxies detected in 3 bands ( 
Lir > 12.22 ). The four infrared sed model compo- 
nents used in Figs 26-29 are the same as those used 
in count models by Rowan- Robinson (2001). 

For all ELAIS galaxies detected in two or more 
bands and with spectroscopic or photometric redshifts 
(a total of 306 galaxies) we have selected the best- 
fitting of the four model components cirrus, M82, 
Arp220 and AGN dust torus, and estimated the bolo- 
metric infrared luminosity (from 3-1000 /xm: here L 
denotes the logio of the luminosity). Wo have used 
the 15 /um fiux (or upper limit) to estimate the lumi- 
nosity in an AGN dust torus component, Ltor, if the 
15 fj,m emission is interpreted as due to such a com- 
ponent (without corroboration at longer wavelengths, 
Ltor gives a more conservative estimate of luminosity 
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Figure 19. 15/^m/20cm colour versus z for ELAIS Cata- 
logue sources, with loci for cirrus (C, solid line), starburst 
(B, dotted line) and Arp220 (A, broken line) seds. Filled 
circles: AGN, open circles: galaxies. 



Figure 20. 15/6.7 fim, colour versus z for N2 and SI, 
with loci for cirrus (C, solid line), starburst (B, dotted 
line), Arp220 (A, short broken line) and AGN dust tori 
(Q, long broken line) seds. 



than starburst or cirrus models). We estimated the 
corresponding optical bolometric luminosities (0.04-3 
^m), using the photometric sed templates (Rowan- 
Robinson 2003). Figure 30 shows a plot of {Lir ~ Lopt) 
versus Lir. For ease of comparison with earlier work, 
an Einstein de Sitter model (fi = 1, A = 0) model, 
with Ho = 50, is used. 

Galaxies with lower Lir and Lir — Lopt < 0.0 are 
prcdomincutly fitted with cirrus models, but cirrus 
models are also required at higher luminosities and 
values of Lir — Lopt as high as 1.0. This would seem 
to be consistent with the concept discussed by Efs- 
tathiou and Rowan-Robinson (2003) of cirrus galaxies 
with 1 < Av < 3 and high total luminosities, as an 
alternative explanation of many of the galaxies de- 
tected in submm surveys. It is also consistent with 
the source-count models of Rowan- Robinson (2001), 
in which the quiescent, cirrus galaxy component un- 
dergoes the same strong luminosity evolution as star- 
burst galaxies. This high incidence of cold, very lu- 
minous, galaxies with high infrared-to-optical ratios 
is one of the main new results of the ELAIS survey. 
The expected trend towards M82 and Arp220 -type 
infrared seds as we go to higher infrared luminosities 
and higher values of Lir — Lopt is also clearly seen. 

Figure 31 shows a similar plot for galaxies for 
which an AGN dust torus component is preferred. 
Where gala^xies are detected only at 6.7 and 15 
/xm, and starburst luminosities in the hyperluminous 
range are implied, we have always prefered a dust 
torus model as the more conservative assumption. For 
Type 1 AGN, L be interpreted as the log 




2.5 



Figure 21. 90 //m/r-band colour versus z for ELAIS 
galaxies, with loci for cirrus (C, solid line), starburst (B, 
dotted line) and Arp220 (A, broken line) seds. 



of the covering factor. The observed values are in the 
range -0.7 to 0, implying covering factors in the range 
20-100 %, with a median value of 30 %. Values greater 
than would imply that the optical QSO must have 
suffered some extinction. 
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Figure 30. Lj^ — Lopt versus Lj^ for ELAIS galaxies detected in 2 or more ISO bands and with spectroscopic or photo- 
metric redshifts. Filled circles: cirrus sed, open squares: M82 starburst, open triangles: Arp 220 starbursts, crosses: M82 
stajrburst+AGN dust torus. 



7 ULTRALUMINOUS AND 

HYPERLUMINOUS INFRARED 
GALAXIES 

We have found a surprisingly large population of 
Arp220-like ultraluminous infrared galaxies {Lir > 
12.22) in this survey (14 % of the 15 /xm galaxies 
with known z). Morel ot al (2001) reported the first 
hyperluminous infrared galaxy in the survey. In fact 
there are a total of 9 candidate hyperluminous galax- 
ies {Lir > 13.22, Rowan-Robinson 2000) in the sur- 
vey, listed in Table 4 (note that for galax;ies ir sed- 
typc — 6, the luminosity in the starburst and AGN 
dust torus components must be added to get the to- 
tal infrared luminosity). A further 45 galaxies would 
be hyperluminous if their 15 irni emission were inter- 
preted as due to a starburst sed (rather than the more 
conservative assumption that it is due to AGN dust 
torus emission). All 9 objects in Table 4 appear to 
be quasars: selection at 15 /xm does favour detection 
of AGN dust tori. For 4 galaxies this classification is 
based only on a photometric redshift, and this needs 
to be confirmed by spectroscopy (especially where Zph 



= 2-2.5, since there is a strong possibility of aliassing 
- see section 5). The large number of ultraluminous 
and hyperluminous infrared galaxies is probably a re- 
flection of the very steep rise in the star-formation 
rate between z = and z = 2, and the associated 
strong evolution in the AGN population (cf Pozzi et 
al 2003, LaFranca et al 2003), though gravitational 
lensing could also play a part. 

We also find 9 ERO's in the survey, defined as (r- 
K) > 6, all from the K-band photometry of Vaisanen 
et al (2002) and Rigopoulou et al (2003, in preparar 
tion), and these are listed in Table 5. 3 of the objects 
have photometric redshifts and have the scds of el- 
liptical galaxies at a redshift ~ 1. Predicted (r-K) 
colours as a function of z show that elliptical galaxies 
have values > 6 for 1 < x < 4, 

In the optical and near infrared the seds of dusty 
starbursts like M82 and Arp 220 look very like ellip- 
ticals due to the extinction of the young stellar com- 
ponent. Table 5 gives the values of logio{Si5 / Sr) and 
logioiSrad/Sr), which are consistent with all 9 objects 
being highly extinguished starbursts like Arp 220, at 
« ~ 1 (cf Arp 220 model curves in Figs 16 and 17). 
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Figure 22. 175/^m/r-band colour versus z for galaxies in 
Nl and N2, with loci for cirrus (C, solid line), starburst 
(B, dotted line) and Arp220 (A, broken line) seds. 

However only a small fraction of EROs will fall into 
this category of highly extinguished starbursts. None 
of the 17 EROs with (R-K) > 6 in Table 2 of Roche 
et al (2002), which are located within a small 81.5 
sq arcmin area of N2, are detected by ISO (one is 
an ELAIS radio source, 163657+410021). Counts of 
EROs in the Nl and N2 areas, and estimates of their 
space-density, are given by Vaisanen and Johansson 
(2003). They conclude that the redshifts of their ERO 
sample lie in the range 0-7-1.5 and estimate the frac- 
tion of strong starbursts ( > iQMo/yr) to be < 10%. 



8 DISCUSSION AND CONCLUSIONS 

We have presented the Final ELAIS Catalogue at 
U,g',r',i',Z, J,H,K, 6.7, 15, 90 and 175 /xm, and 20 cm. 
The process of band-merging and optical association 
of the sources has given considerable insight into the 
cxtragalactic infrared populations present. Although 
the largest single population is relatively nearby (z < 
0.2) cirrus galaxies, there is a surprisingly large pop- 
ulation of ultraluminous infrared galaxies (14 % of 
15 /jm galaxies), many of them highly obscured star- 
bursts like Arp 220. There is also a significant pop- 
ulation of luminous (Lj^ > 11.5) cool galaxies, con- 
sistent with the idea that the quiescent component of 
star-formation in galaxies undergoes the same strong 
evolution &a the starburst component. There appear 
to be a small proportion of genuinely optically blank 
fields to r' = 24 (10 % at 15 /im, 3 % at 6.7 /xm, up to 
20 % at 90 /ttm, and 1 % at 175 fim) which must have 
high infrared-to-optical ratios, be at higher redshift (z 



Figure 23. 175/90 /tm, colour versus z for galaxies in Nl 
and N2, with loci for cirrus (C, solid line), starburst (B, 
dotted line) and Arp220 (A, broken line) seds. Horizon- 
tal bars at right indicate dust colour temperatures from a 
uB^{T) model. 

> 0.6), and have high infrared luminosities. They are 
therefore dusty luminous starbursts or Type 2 AGN. 

The ELAIS survey provides a strong basis for 
the SIRTF surveys. The SWIRE Legacy Survey will 
include the ELAIS Nl, N2, and SI areas. It will be 
worthwhile to obtain better seds and ir spectroscopy 
for the luminous infrared galaxies in the survey. As 
SIRTF does not have a 15 /xm band, the ELAIS data 
will be a useful complement to the SWIRE data. The 
results presented hero illustrate the importance of 
multiband optical data for photometric redshift de- 
termination. The far infrared colour-colour diagrams 
discussed here illustrate the difficulties of trying to 
determine redshifts from far infrared data alone. 
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Table 4. Hyperluminous infrared galaxies in the ELAIS Catalogue 



name 




fwFS 


Tltyp Zph 


^spect 






Lopt 


"ISO 






Ltor 


ELAISC15-J002925-434917 


18.63 






3.094 


5 


3 


13.59 


1 


2 


14.80 


13.54 


ELAISC15-J003213-434553 


17.09 






1.707 


5 


3 


13.57 


1 


2 


13.97 


13.29 


ELAISC15-J050152-303519 


17.86 




4 0.0233 


1.813 


5 


4 


13.62 


2 


4 


13.25 


13.07 


ELAISC15-J160419.0+541524 


18.00 


0.99 


8 2.548 








13.42 


1 


2 


14.51 


13.39 


ELAISR 160758+542353 


22.70 


0.10 


8 2.311 








11.48 


2 


2 


13.80 


-12.57 


ELAISC15-J161259.2+541505 


19.06 


0.98 


8 2.548 








12.98 


2 


6 


13.76 


13.03 


ELA1SC15-J163739. 2+405643 








1.438 


5 


1 




3 


2 


13.42 


12.67 


ELAISC15-J164010.1+410521 


16.95 


1.00 


8 1.399 


1.0990 


5 


6 


13.17 


3 


6 


13.39 


12.92 


ELAISC15-J164018.4+405812 


18.06 


0.99 


8 2.311 








13.33 


1 


2 


14.33 


13.30 



Table 5. EROs in the ELAIS Catalogue 



name 


r' 


K 




^phot 


IgiSis/Sr) 


IgiS 


ELA1SR160721+544757 


24.11 


18.076 








3.60 


ELA1SC15-J160913. 2+542320 


24.29 


17.180 






3.48 




ELA1SR161030+540247 


23.53 


16.400 








2.91 


ELAISR161046+542329 


23.15 


17.090 


2 


1.138 




2.28 


ELAISC15-J163536.6+404754 


23.38 


17.340 


1 


0.995 


2.71 




ELAISR163555+412233 


23.29 


16.789 


1 


1.291 


2.62 




ELAISR163723+410526 


23.71 


17.449 








3.33 


ELAISC15-J163748.1+412100 


24.50 


17.926 






3.14 




ELAISR163758+41 1741 


23.79 


17.413 








2.32 




Figure 24. U — r' versus versus z for galaxies in Nl and 
N2. Model loci arc for an E galaxy (dotted line), Sbc 
galaxy (long-dashed line) and starburst (dashed line). 

Blommacrt ct al, 1998, ISOCAM Handbook 

Chapman S.C., Small L, Ivison R.J., Helou G., Dale D.A., 
Lagache G., 2002, ApJ 573, 66 

Cihegi R, McMahon R.G., Miley G., Gruppioini C, 
Rowan- Robinson M., Cesarsky C, Danese L., Prances- 
chini A., Gcnzcl R., Lawrence A., Lemke D., Oliver S., 
Puget J.-L., Rocca-Volmerange B., 1999, MN 302, 222 

Cram L., Hopkins A., Mobasher B., Rowan-Robinson M., 



Figure 25. U — r' colour versus z for star-like objects in 
Nl and N2. Filled triangles denote objects with AGN scds, 
Filled circles arc spectroscopic AGN (Sy 1 and 2) and open 
circles are objects with galaxy seds. The model loci arc for 
an AGN (solid line), E galaxy (dotted line), Sbc galaxy 
(long-dashed line) and starburst (dashed line). 

1998, ApJ 507, 155 
Crampton D., Cowley A.P., Hartwick F.D.A., Ko P.W. 

1992 Astron. J. 104, 1706 
Dole H., Gispert R., Lagache G., Puget J.-L., Bouchet 

F.R., Cesarsky C, Ciliegi P., Clements D.L., Dennefeld 
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Figure 26. 9 ELAIS galaxies detected in all 4 ISO bands 
for which a standard cirrus model provides an excel- 
lent fit to the spectral energy distribution. Filled cir- 
cles: ISO data, filled triangles: IRAS data. From bot- 
tom: ELAISC15 163359.2-1-405303 (Li^=10.13, ntyp=l), 
163401.8-1-412052 (10.42, 3), 163506.1-1-411038 (11.01, 
(4)), 163525. H-405542 (10.46, 4), 163546.9-1-403903 
(11.19, 3), 163548.0-1-412819 (11.21, 5), 163607.7-1-405546 
(11.41, 4), 163613.6-1-404230 (10.67, 1), 163641.1+413131 
(11.04, 2). 
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Figure 27. 9 luminous ELAIS galaxies detected in 3 
ISO bands for which a standard cirrus model provides 
an excellent fit to the spectral energy distribution. Filled 
circles: ISO data, filled triangles: IRAS data. From bot- 
tom: ELAISC15 050225-3041112 (Li^=11.59, ntyp=l), 
160443.5-f543332 (12.30, 1), 160552.5-1-540650 (11.99, 3), 
160553.3-1-542225 (11.52, 1), 160945.7-1-534944 (11.68, 1), 
161041.2-1-541032 (11.51, 4), 163242.4-1-410847 (12.11, 2), 
163449.5+412048 (11.61, 1), 163741.3+4111913 (11.97, 
3). An additional star burst component is shown for 
160552.5+540650 (dotted curve). 
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Figure 28. 6 ELAIS galaxies detected in 3 or 4 ISO 
bands for which an Arp200 model provides an excel- 
lent fit to the spectral energy distribution. Prom bottom: 
163412.0+405652 (11.40, 2), 163608.1+410507 (12.16, 1), 
160934.7+53220 (11.46, 1), 161300.8+544838 (11.98, 2) , 
163449.5+412048 (12.15, 1), 163708.1+412856 (11.77, 1). 
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Figure 29. 5 ELAIS galaxies detected in 3 ISO bands 
compared with an M82 starburst model. For two of the 
galaxies an AGN dust torus component has been in- 
cluded to fit the mid-ir data (model from Rowan- Robinson 
1995). The top two loci are sources identified with quasars, 
with evidence only for dust torus emission in the in- 
frared. From bottom: ELAIS 163751.4+413027 (12.20, 4) , 
164010.1+410521 (13.17, 8), 160443.5+543332 (12.41, 1), 
163615.7+404759 (13.00, 2), 163741.3+411913 (12.11, 3), 
163352.4+402112 (11.79, 7), 163502.7+412953 (11.48, 7). 




Figure 31. Li^—Lopt versus Lir for ELAIS AGN dust tori 
detected in 2 or more ISO bands and with spectroscopic 
or photometric redshifts. Filled circles: Sy 1 (from optical 
spectroscopy) , open circles: Sy 2, or galaxy sed in optical. 
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